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Abstract. We explore possible constraints on the inflationary equation state: p = wp. While 
w must be close to —1 for those modes that contribute to the observed power spectrum, for 
^ those modes currently out of experimental reach, the constraints on w are much weaker, 

~~* with only w < —1/3 as an a priori requirement. We find, however, that limits on the 

reheat temperature and the inflationary energy scale constrain w further, though there is 
^f) still ample parameter space for a vastly different (accelerating) equation of state between 

[^ the end of quasi-de Sitter inflation and the beginning of the radiation-dominated era. In the 

^P event that such an epoch of acceleration could be observed, we review the consequences for 

the primordial power spectrum. 
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1 Introduction 

The theory of inflation has the twin properties of explaining away (at least to some degree) 
many of the puzzles of the very early universe in a broad class of general models, whilst 
simultaneously remaining remarkably hard to pin down in specific detail^. The key success 
of the inflationary theory is its prediction of an almost scale invariant power spectrum of 
primordial fluctuations [2-6], a prediction that is borne out by what we measure of the 
perturbations in the cosmological fluid today. However, the limited information we can 
extract about this power spectrum means that a wide class of inflationary models are viable. 
Moreover, the particular nature of that class of models and the constraints upon them are 
predicated on many assumptions about both the nature of inflation and the evolution history 
of the universe. 

Whilst inflation is typically modeled by a (possibly multi-component) scalar fleld, the 
essential feature is that the expansion of universe undergoes a temporary accelerating period, 
where a > 0. During accelerated expansion fluctuations "exit" the horizon (i.e. become larger 
than the physical size of the temporary event horizon) . Then, once the expansion begins to 
decelerate, these fluctuations "re-enter" the horizon in reverse order from their exit. This 
means that the smallest scale perturbations we observe were produced last, closest to the 
start of decelerated expansion. 

The top graphic in figure 1 schematically illustrates this, with a logarithmic plot of the 
quantity k^,/aH. This quantity characterizes the size of the comoving horizon relative to a 
particular scale, /c=k, with k^,/aH = 1 corresponding to the scale exiting (or re-entering) the 
horizon. The slope of the line is related to the dynamics (or, in the language we will be 
using presently, the equation of state) of the cosmological fluid - when the slope is nega- 
tive(positive), the expansion is accelerating(decelerating). 

In this paper we eschew the aforementioned scalar field paradigm, and rather consider 
the equation of state, bounding, as best we can, the equation of state parameter w, defined by 
p = wp. The Friedman equations (along with the null energy condition) mean that a positive 
a implies that —1 < w < —1/3. Changing the equation state changes the tilt of of the slope of 
the line in figure 1 . We are primarily interested in changes to the negative part of the slope 
(i.e. the accelerated expansion) that are consistant the current observational constraints. 
One could also consider bounds on the equation of state during decelerated expansion, which 

^For an alternative solution to some of these puzzles see [1] 






Figure 1. Schematic evolution of the quantity k/aH for a particular scale, /c*. Changing the equation 
of state before(after) the end of inflation changes the negative(positive) part of the slope 



has been done by, for example, [7] and [8]. The two possibilities are schematically shown 
respectively on the left and right of the lower half of figure 1. As we shall see, the broad cause 
of the uncertainty in the inflationary equation of state (and any uncertainity in the equation 
of state immediately after inflation) it that we only have a direct probe of the primordial 
spectrum for a subset of the e- foldings that correspond to our current observable universe. 

An equation of state analysis of inflation has been carried out recently by Ilic, Kunz, 
Liddle and Frieman in [9], our approach is somewhat different, in particular we explore 
possible constraints on w outside of those directly imposed by the power spectrum. Said 
power spectrum is almost scale-invariant from the size of the observable universe down to 
scales of around a Mpc; as the authors show in [9], this places tight bounds on the deviation 
of w from —1. However, as we will describe below, w is less tightly constrained for those 
periods of accelerated expansion corresponding to shorter length scales. 

We explore the possibility of a significantly varying w before the start of decelerating 
expansion by dividing the inflationary epoch into two distinct periods. The first, which for 
clarity we will refer to as "inflation", has tu ~ — 1 and is responsible for producing the nearly 
scale-invariant power spectrum we see on large scales. The second, from now on we call this 
"accelerated expansion" , lasts as long as is necessary to be consistent with particular choices 
of the reheat scale (where deceleration begins at reheating) and the inflation-acceleration 
transition scale. 

In section 2 we review the derivation of the mode functions, and find an expression for 
Us in terms oiw. In section 3, we outline our model of the thermal history of the universe and 
discuss how constraints on observable wavenumbers in the power spectrum serve to bound w 
(and thus Ug) during a period of accelerated expansion. We finish with some brief concluding 



remarks on what can be learned from this approach. 

2 Fluctuations During Accelerated Expansion 

2.1 Not Power-Law Inflation 

Before marching onwards with our calculations, we should note that periods of expansion 
with w 7^ —1 (but < —1/3) give rise to what is more commonly known as power-law inflation, 
where the scale factor grows like t^ {p > 1 is a function of w ~ the detailed solution can be 
found below). This is a well studied situation [10-12] in inflationary cosmology; it arises, 
for example, when we have a single scalar field with an exponential potential. We should 
note, however, that our setup is different in two important ways. Firstly, we are positing 
accelerating expansion in addition to inflation - i.e. it is inflation that provides the nearly 
scale-invariant power spectrum we need to explain CMB and matter power spectrum data. 
All we allow our period of accelerating expansion to do is to provide additional e-foldings 
as needed to solve the horizon problem and to change (possibly) the power spectrum at the 
highest wavenumbers. 

Furthermore, our treatment of the quantum vacuum for the modes produced during 
accelerated expansion is different. Rather than assuming a Bunch-Davies vacuum state [13, 
14], we instead take the modes during inflation as an initial condition on which to match 
those produced during accelerated expansion - this is equivalent to assuming a sudden change 
in the equation of state from w ~ —1 to —1 < w < —1/3. For obvious reasons (i.e. we match 
only the shortest wavelength modes) this gives the same answer as assuming that the Bunch- 
Davies vacuum is the appropriate one to use for power law accelerating expansion. In the 
case of power-law inflation for sufficiently large powers {w sufficiently close to —1) one is 
close enough to de Sitter space such that the de Sitter invariant a- vacua give an appropriate 
family of privileged states from which the Bunch-Davies vacuum is chosen. For larger values 
of w (smaller powers), while one can construct the Bunch-Davies state, it no longer belongs 
to such a privileged invariant set, and thus is no longer an obvious choice for an appropriate 
vacuum. 

Caveats and qualifications out the way, let us now move on to calculating the fluctuation 
spectrum. 

2.2 Mode Functions 

We are interested in the scalar fluctuations produced during an epoch of accelerated expansion 
(with the equation of state is given hy p = wp and —l<w< —1/3). Following standard 
treatments [15, 16] we have the gauge invariant quantity, IZ, which for a single scalar field is 
conserved outside the horizon during inflation {IZ = ^ + H6(j)/(f)). The Mukhanov variable 
is: 

v = zn, (2.1) 

with z is defined in terms of the background scalar field as: 

(2.2) 
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In terms of the equation of state parameter, z^ = ^^a'^i^ + w)- The Fourier transform of v 
satisfies the Mukhanov equation [6, 17]: 

(Pvu /,9 ld?z\ 

+ {k^--:r^]^k = o, (2.3) 
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where r is the conformal time. 

Solving this equation allows us to find the mode functions during both inflation and our 
second period of accelerated expansion. This, in turn, will let us find an expression for the 
spectral index. We summarize the results below. 

Solutions 

During inflation, with w = Wmi ~ —1 and a = afS f^*'~^f' , the mode function is given by (in 
the Bunch-Davies vacuum): 
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From equation 2.1, this give IZt as: 
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During an epoch of accelerated expansion (matching the scale factor and energy density 
to the prior inflationary period) the scale factor evolves as: 
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with H given by: 



H = Hf 



a/\f(i+™) 



and T defined in terms of a by: 
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This allows us to write equation 2.3 as: 
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where Tf = — ^ ^ . The general form of the solution is given by a linear combination of 
Hankel functions: 
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(2.10) 



(2.11) 



Then 7^ during this period is given by: 
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The functions A{k) and B{k) can be found by matching this solution with that given 
in equation 2.5 - this, of course, assumes an instantaneous transition from inflation to a 
second period of acceleration. Doing this for all values of k is non-trivial, in particular for 
k corresponding to the horizon size, only numerical evaluation of A{k) and B{k) is possible. 
However, for modes asymptotic in \kTf\ it is relatively straightforward to find the curvature 
perturbation. Far outside the horizon, \kTf\ « 1, the curvature perturbation is "frozen-in" 
and curvature perturbation remains at the constant value set by the inflationary epoch. Deep 
inside the horizon, \kTf\ >> 1 matching solutions with equation 2.5, we obtain (for r > tj): 
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(2.13) 
While it is relatively easy to extract observables (in particular the value of Ug) from the 
above expression (we do so below), we should note that if a secondary period on inflation 
produced visible consequences, the largest scales at which such effects could be observed 
do not correspond to either of the two limits discussed above. Rather, since inflationary 
perturbations would be responsible for large scale perturbations that we already measure, 
from the size of the observable universe down to scales less than ~ IMpc, the largest scales 
from a secondary period of accelerated expansion would be those produced right after the end 
of inflation, when krj ~ 1. At these scales, we have to match our solutions numerically. We 
have carried out numerical investigation of the matching for krf ~ 1 , and flnd that the form 
of Vk approaches the asymptotic results within a few percent for 0.5 < krj < 2. Accordingly, 
we concentrate on the form of the mode function given in equation 2.13 in order to obtain an 
analytic expression for the primordial power spectrum of a secondary period of acceleration, 
though we note that a more precise numerical calculation would be needed in the event of 
an actual measurement. 

The Spectral Index rig 

In order to construct a physical observable, we calculate the scalar spectral index. To do 
this we take the superhorizon limit, \kx/Q\ << 1,^ and use the asymptotic behavior of the 
Hankel function for \y\ « 1 (this result, of course, only applies to modes which also have 
\kTf\ >> 1 - i.e. modes which were well inside the horizon at the end of inflation): 

fl(..^)[,|^^!lM(|)-". (2.14) 



From this we obtain: 
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Therefore, |7^fcp oc k ^", and the scalar spectral index, n^ (deflned through \TZk\'^ oc k '^+"''), 
is given by: 

ns - 1 = 3 - 2ra = / '^ { . (2.16) 



^One can check this is indeed the superhorizon hmit by finding expressions for both the physical wavenum- 
ber, k/a, and the event horizon, a j^ ^,2^^, ,. da', in terms of the quantities Q and x defined in equation 
2.11. 



Clearly, for —1 < w < —1/3, we have 1 > Ug > — oo. Having found the relationship between 
w and Us our next goal is to see how constraints can be brought to bear on these quantities. 
We do this by considering the thermal history of the universe. 



3 Constraints Prom Thermal History 

As discussed above, fluctuations produced during a period of accelerated expansion exit 
the horizon, and when they do so the associated curvature perturbation does not evolve 
(assuming that a single component fluid is driving the acceleration), until it re-enters the 
horizon at some later time, when the expansion is decelerating. Since our posited period of 
additional accelerated expansion takes place after inflation, it will affect the power spectrum 
at shorter scales. By considering the evolution history of the universe and the various bounds 
observations of the power spectrum place on said history, we can constrain some of the 
parameters of any additional period of accelerated expansion. 

We assume a sequence of epochs after inflation, with a variety of different equations 
of state: during inflation we have w = —1, during our accelerated period of expansion we 
have —1 < Wa < —1/3, then from t = t^h (the time of reheating) the universe evolves 
in the conventional manner as an appropriate mixture of radiation, matter and A (this, of 
course, assumes an instantaneous reheat). Solving the Friedmann equations in the various 
different epochs and matching the scale factor and energy density across boundaries gives the 
solutions shown in table 1 (for clarity we've glossed over the details of scale factor evolution 
after reheating, but this will be taken into account below when finding constraints). 
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Table 1. Cosmological solutions for epochs of interest. These equations are somewhat over- 
parameterized, all parameters can be expressed in terms of pf and various scale factor normalizations. 

We can use the above history to give us a bound on the combined number of e- 
foldings from inflation and our accelerated period of expansion. Noting that during radiation- 
domination and other decelerating periods of expansion the particle horizon grows faster than 
the universe (and thus faster than the perturbations in the cosmological fluid), we see that 
in the far past modes that are currently inside the horizon were far outside it. This is the 
horizon problem, and it is solved by having a period of expansion where the horizon grows 
slower than the universe, which lasts long enough to ensure that scales within horizon today 
were also within the horizon in the early universe. The number of e-foldings that are required 
to solve the horizon problem is fixed by the amount of horizon growth relative to the size of 



the universe since deceleration began [18], or: 
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h is the Hubble constant in units of 100km s~^ Mpc^ . Furthermore, the number of e-foldings 
that took place between the time at which the mode corresponding to the size of the horizon 
today previously held that honor, and the time of reheating is given by: 
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(3.2) 



tqQ is the time at which the mode corresponding to the size of the horizon today was produced. 
Since for our model this mode was produced during inflation, and since H is approximately 
constant during the inflationary epoch, we can replace tq^ with tf. 
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Na is the number of e-foldings from accelerated expansion. Thus we have: 
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(3.3) 



(3.4) 



Incidentally, since Nobs ^ ^o and since the observed scale-invariant spectrum of per- 
turbations is at scales from the current size of the observable universe down, requiring that 
the primordial power spectrum (at large scales) is explained by inflation guarantees that 
the horizon problem will be solved. It's clear from equation 3.1 that the minimum required 
number of e-foldings is fewer for lower scale reheating. We cannot, however, push this scale 
arbitrarily low - we know, for example that the universe was radiation dominated at the time 
of big bang nucleosynthesis, which corresponds to a scale of around 1 MeV (in [20] Hannes- 
tad finds T^h > 4 MeV or > 1 MeV if reheating direct to neutrinos). We thus work with a 



conservative lower bound on the reheat scale of p, 
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10 MeV or - this would give > 19 



visible e-foldings (the precise bound depends also on the value of Ni)'^ . This analysis, and 
that which follows below, assumes instantaneous reheating - if this were not the case the true 
value of Na would be somewhat smaller, with the rest of parameters adjusted accordingly. 

In our scenario the e-foldings within the observable universe are divided between in- 
flation and accelerated expansion. We can place a lower limit on the portion of visible 
e-foldings coming from inflation by considering the form of the power spectrum, in particular 
the size of the deviation from scale invariance. To do this we use a recent reconstruc- 
tion of the power spectrum from Peiris and Verde [21]. The authors reconstruct ns{q) for 
0.0001 < q [/i/Mpc] < 3, finding that 0.7 < ns{q) < 1.3 - the values of this bound depend on 
the details of the reconstruction, and it is somewhat tighter for q away from the upper and 
lower limits. Let us define qmax as the physical wavenumber corresponding to the smallest 
scale perturbation produced by infiation and go as the physical wavenumber corresponding 
to the horizon today, then the number of visible e-foldings produced by inflation is given by: 



iV,. = ln 



(3.5) 



One might complain that baryogenesis places tighter bounds on Trh\ while this is true for most models of 
baryogenesis, it can be avoided in the AfHeck-Dine scenario [19]. 
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Figure 2. The left-hand plot shows the relationship between Na and w for different choices of pj^^ 
(with TV, - In W): 10^^ MeV for the uppermost dotted hne, lO^^ j^gy for the dashed line, 10* MeV 
for the dot-dash line, 10'' MeV for the widely spaced dashed line and 10 MeV for the solid line. The 



right-hand plot shows the same for differing Ni (with p. 



1/4 
rh 



10* MeV): In 10'' for the uppermost 



dotted line, In 10^ for the dashed line, In 10^ for the dot-dash line. In 10^ for the widely spaced dashed 
line and In 10* for the solid line. 



Then from equation 3.4 and noting Nobs = Na + Ni, we have: 
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(3.6) 



Na can also be expressed in terms of the ratio of the scale factors at the beginning and end 
of accelerated expansion: 
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Bounding pV^ < 10^^ GeV (this bound comes from the amplitude of CMB fluctuations, 10~^ 
and the tilt-imposed upper limit of the slow roll parameter, e < 0.05, see e.g. [18] for details) 
we can use equations 3.6 and 3.7 to find allowed values of Na and Wa for different choices of 

Prh- 

In figure 2 we illustrate the nature of the Na-Wa relationship as we vary pj^^ (left-hand 
plot) and Ni (right-hand plot); note that pf varies as one moves along the curves, each of 
which terminates when the upper limit for pf is reached. It is clear that for a sufficiently 
low reheat temperature our current observations of the power spectrum could be consistent 
with having a period of quasi de Sitter inflation followed by a period of accelerated expansion 
with an equation of state significantly different from w ~ —1. Moreover, such a situation 
can be arranged with Ni sufficiently large such that the perturbations sourced by accelerated 
expansion are at length scales much smaller than those currently probed by experiment. On 



the other hand, if we wish to push up as close as possible to the current measurement of 
Ni ~ In 10^ and thus be on the boundary of having an observable signature of accelerated 
expansion, we find the most dramatic result we can have (i.e. the largest w), comes from 
having a minimal reheat temperature, and would give us ~ 93 e-foldings of accelerated 
expansion with w ~ —0.41 and Ug ~ —14.7. 

As our ability to measure the power spectrum on the shortest possible length scales 
improves, it is possible that we may be able to see the effects of a variation in the equation 
of state. However, given the difficulty of extracting information about the primordial power 
spectrum from scales much smaller than the scale at which local gravitational effects (i.e. not 
cosmology) dominates the physics it seems unlikely that any distinct signature of a second 
period of acceleration, unless its effects start at scales close to those at which we are currently 
confident of a scale-invariant spectrum. This restriction applies not only to the tilt, but also 
to other potential probes of the equation of state that rely on direct knowledge of the the 
spectrum (e.g. the running of Ug, higher order correlations etc.). 

Although direct measurement of the spectrum is progressively harder at shorter and 
shorter distance scales, there are indirect probes that could be useful bounding the amplitude 
of primordial fluctuations. In particular, we know there must have been sufficient power at 
short distances to give rise to star production at an early enough times to produce the 
reionizened universe that we observe for z < 6 (some of the relevant issues are discussed in 
[22]). In principle this requirement may conflict with the most extreme scenarios we consider 
above, since a strongly red-tilted spectrum means that considerably less power is present 
at short scales when compared to the 6p/p ~ 10~^ that we have at large scales. However, 
although such constraints depend on the details of reionization in a manner beyond the scope 
of this paper, they seem unlikely to provide much stronger bounds on Na and Wa for any 
given qmax/qo and p^^ . 

Before moving on to our conclusions, let's consider what might happen if there was a 
secondary period of acceleration with a more complicated equation of state than that given 
by a constant w. To understand how such a period (distinguishable for inflation in principle) 
can be constrained in general it is helpful to consider again the relationships Na = N^i^g — Ni 
and Nobs > -^o (the latter is implied by equations 3.3). Together these give: 
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Given the lower bounds on the reheat temperature and the number of e-foldings of almost 
scale-invariance it's clearly possible to have at least a few e-foldings of arising from a second 
accelerating epoch (with 0(10) MeV reheating and scale invariance down to a Mpc, there 
are 8 unobserved e-foldings). Higher reheating scales would leave more possible e-foldings 
hidden. The precise nature of the accelerating expansion that is responsible for those e- 
foldings outside our current observations is almost free of constraints, subject only to the 
requirement that the energy density must change by a sufficient amount between the end of 
inflation and the start of decelerated expansion after reheating. This leaves a large set of 
possibilities available for an alternative acceleration mechanism between the end of inflation 
and the start radiation-domination, though one would expect that unless the equation of state 
parameter is changing rapidly compared to the Hubble constant, the space of possibilities 



would be constrained in a similar manner to the constant w case discussed here^. 

4 Conclusions 

In this paper we have explored the possibilities of a secondary period of accelerated expansion 
in the universe's history. With such a period immediately following inflation we show that the 
bounds imposed by considering current measurements of the power spectrum are relatively 
weak - with reasonable reheat temperatures and inflationary scales it is relatively easy to 
fit in a secondary period of accelerated expansion with a vastly different equation of state 
to inflation and still have a nearly scale-invariant spectrum at scales from the size of the 
universe down to a Mpc and below. 
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